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Abstract 

There may be non-thcmial cosmic rays during big-bang nucleosyn- 
thesis (BBN) epoch (dubbed as BBNCRs). This paper investigated 
whether such BBNCRs can be the origin of Lithium problem or not. 
It can be expected that BBNCRs flux will be small in order to keep 
the success of standard BBN (SBBN). With favorable assumptions on 
the BBNCR spectrum between 0.09 - 4 MeV, our numerical calcula- 
tion showed that extra contributions from BBNCRs can account for 
the ^Li abundance successfully. However ^Li abundance is only lifted 
an order of magnitude, which is still much lower than the observed 
value. As the deuteron abundance is very sensitive to the spectrum 
choice of BBNCRs, the allowed parameter space for the spectrum is 
strictly constrained. We should emphasize that the acceleration mech- 
anism for BBNCRs in the early universe is still an open question. For 
example, strong turbulent magnetic field is probably the solution to 
the problem. Whether such a mechanism can provide the required 
spectrum deserves further studies. 
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1 Introduction 



Big Bang cosmology is an excellent model to describe our Universe [1]. Al- 
beit its success there is still puzzling "Lithium problem" [2, 3] in its building 
block - big bang nucleosynthesis (BBN). Namely the theoretical predicted 
primordial ^Li abundance is higher than observation while the ^Li abundance 
is lower. Primordial ^Li abundance from measurements of metal-poor halo 
stars is ^Li/H= (1 ~ 2) x 10"^° [2]-[9], while prediction by the standard BBN 
(SBBN) model is three to five times higher: ^Li/H= (5.24l°;^^) x 10-^° [10]. 
Meanwhile ®Li abundance from observation is ^Li/B.^ 6 x 10~^^ [3], a factor 
of about 1000 higher than the SBBN model prediction [11]. 

Extensive investigations have focused on Lithium problem. The discrep- 
ancy could be due to astrophysical origin [12, 13, 14]. While it is also possible 
that the discrepancy is arising from Physics beyond the SBBN model. Inves- 
tigations showed that Lithium abundance would change by varying the nu- 
cleon effective couplings as well as the mass parameters, namely the neutron 
lifetime, the neutron-proton mass difference and the deuteron binding energy 
etc. [15, 16, 17]. These parameters can be modified due to virtual effects 
of Physics beyond the Standard Model of particle physics. Another possi- 
bility of Physics beyond the SBBN could contain new particles/resonances. 
Such new particles/resonances participate in thermal nuclear reactions during 
BBN [18]-[28], as such Lithium abundance will depend on the detail assump- 
tions about the new particles/resonances. Besides these two approaches, non- 
thermal electromagnetic [29] and hadronic [30]- [35] energy injection was also 
investigated. Non-thermal particles originate from the nuclear reactions be- 
fore the final state particles are thermalized. The preliminary studies showed 
that the mono-energetic injection is hardly the origin of Lithium problem. In 
this paper we will extend, in some sense, such idea to include the non-thermal 
energy injection from cosmic rays during and/or shortly after the epoch of 
BBN. 

In this paper, we propose a possible solution to the Lithium problem 
by not going far from the SBBN model, neither due to the astrophysical 
factors nor by introduction of new particles. In order to destroy ^Li ^, we 
include new contributions from non-thermal particles, namely the accelerated 
SBBN particles. As is known that the non-thermal particles accelerated in 

^ Actually we need to destroy ^Be. The relic primordial ^Li mainly (about 90%) comes 
from the decay of ^Be. 
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astrophysical environment today are called cosmic rays, it is expected there 
is such cosmic rays even in the early stage of BBN. We expect that some 
kind of plasma wave or other mechanisms can feed energy to thermal ions 
[36], provided that there is strong enough turbulence before/during BBN 
epoch. As a consequence the particle energy spectrum deviates from thermal 
distribution. Exploring detailed mechanism on how to induce non-thermal 
component is beyond the scope of this paper and which will be the focus 
of the further work. However we propose a toy model to support our basic 
hypothesis in this paper. 

The primordial magnetic fields might be created at some early stage of 
the evolution of the Universe, for example inflation, the electro-weak phase 
transition, quark-hadron phase transition and so on. As investigated by au- 
thors in ref. [37], after electro-weak phase transition the magnetic field build 
up and evolve with the expanding Universe. We can estimate the strength 
of induced electric field through E ^ AB/At x L ^ B H L, where B oc i?"^ 
is the characteristic magnetic field, H is the expansion rate as inverse of the 
characteristic time, and L is the characteristic length of turbulence which is 
given by equations (47) and (49) of [37] (here L = Iq{R/Rq) which includes 
the effect of cosmic expansion). Using the same initial parameters as in ref. 
[37] and extrapolating to O.OlMeV, we get the induced electric field E ^ 
V/m in plasma with temperature O.OlMeV. A charged particle like proton 
undergoing such electric field will gain energy AE ^ rthcr "^jl Q E, where v\\ 
is the particle velocity parallel to the electric field, q is the electric charge of 
the particle, and Tthev{T) is the thermalization time of a high energy (of order 
MeV) nuclear in plasma of temperature T. For O(MeV) protons in plasma 
of temperature O.OlMeV, rthcr(O.Ol) ^ lO'^ s [38], the possible energy gam 
of an accelerated thermal proton is of order 0.1 MeV. During one free time, 
a proton can travel a distance of Lfree ~ 10~^c xlO~^s 3 x lO^cm (c is 
the light velocity) which is just several percent of the characteristic length 
scale. Protons have the opportunity to be accelerated several times and thus 
to gain a couple of MeV energy, though the probability is low. Thus we can 
expect the energy spectrum of high energy particles will not be too hard. 

This proposed toy model needs hydro-magnetic turbulence as the premise. 
We don't know exactly how hydro- magnetic turbulence develops during the 
BBN. However according to ref. [39], helical hydromagnetic turbulence sur- 
vives until 100 eV. We found that the correlation length and magnetic field 
will be larger if we adopt model in ref. [39] instead. Therefore there should 
be a considerable parameter space to support our hypothesis. 
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In this paper we will examine phenomenologically whether such scenario 
can account for Lithium problem. These non-thermal particles are dubbed as 
BBN cosmic rays (BBNCRs). In the SBBN conventional abundance calcu- 
lation, the exothermic reactions are included. BBNCRs contributions to the 
exothermic SBBN reactions will be investigated. Most importantly, we have 
to include the endothermic reactions in order to destroy ^Be by BBNCRs. 
In the SBBN, the endothermic reactions are thought not important. 

This paper proceeds as following. In section 2 we describe the required 
characteristics of BBNCRs, namely the particle species, the energy range and 
energy spectrum and so on. There are many reactions involving BBNCRs 
so we need to select which reactions may be important, and detailed discus- 
sions on these are contained in this section. The necessary formulae for the 
abundance computation by including BBNCRs contributions are depicted in 
section 3. The numerical results are given in this section. Section 4 con- 
tains our conclusions and discussions. Some calculation details are given in 
Appendix A and B. 

2 BBNCRs and Reactions 

In order to estimate the effects of BBNCRs, we need to know the flux, 
the energy range and the energy spectrum of BBNCRs. They are maybe in 
evolution during its work time with the Universe expanding, so the following 
characters can be seen as an average. Obviously we don't have any knowledge 
on BBNCRs, therefore we examine some constraints on BBNCRs. 

Candidate for BBNCRs may be protons, neutrons, nuclei and electrons. 
Neutrons are hard to accelerate for electric neutrality, and electrons are out 
of consideration in SBBN computation, so we focus on nuclei. Because He- 
lium loses energy more quickly through Coulomb scattering [38, 41, 42] due 
to the higher nuclear electric charge [43], we exclude ^He and ^He as the 
components of BBNCRs. We assume that BBNCRs consist of energetic hy- 
drogens, namely protons, deuterons and tritones. In order not to violate the 
success of SBBN, the intensity of each kind of BBNCRs must be much lower 
than that of the corresponding SBBN particles. For simplicity we assume 
that the fraction of BBNCRs is fixed as one single free parameter e. 

Next we examine the energy range of BBNCRs. In order not to change 
deuterium abundance significantly, we must avoid BBNCRs contributions 
from reaction D(p, n)2H. Note that the threshold energy and cross section 
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of this reaction are 3.337 MeV and ~ 10~^ barn respectively. We exam- 
ined the endothermic reactions with threshold energy [44] below 3.337 MeV. 
In order to destroy ''Be by proton, the most effective one is ^Be(p,pa)^He 
with threshold energy Eth =1.814 MeV. Other endothermic reactions besides 
^Be(p,pa)^He are summarized in table 1. 



Table 1: Endothermic reactions 
by energetic hydrogen with threshold energy below 3.337 MeV. 



Process 


Threshold/MeV 


Effect 


^Be((i,3He)6Li 


0.144 


destroy ^Be indeed, but less important 
than ^Be(p,pa)^He, for D is much less 
than p; produce ^Li 


'^U{d,pf Li 


0.247 


destroy '^Li, not important, for we 
need to destroy ^Be 




0.747 


transformation between ^Li and ^Be, 






but less important than ^Li(p, n)^Be 

transformation between T and ^He, and T 

and ^He abundances are assumed to change little 




1.019 


''U{d,t)^U 


1.278 


destroy ''Li; produce ^Li 


^Li{p,pa)D 


1.721 


destroy ^Li, negative, for we 
need to produce ^Li 


^Be(p,pQ;)^He 


1.814 


the most important process to 
destroy '^Be 


^U{t,npYU 


1.850 


destroy ^Li 


'^l-ii{p, n) ''Be 


1.880 


transformation between ^Li 
and ''Be, important 


^U{d,a)2D 


1.967 


destroy ^Li 


^Be{d,daf}ie 

^U(t,da)T 

^Be(t,to)^He 

=^He((i,2p)T 

^Be{d,nfB 


2.041 


destroy ''Be, but less important 


2.213 


destroy ^Li 


2.268 


destroy ''Be, but less important 


2.436 


transformation between T and ^He 


2.688 


destroy ''Be, but less important 


^Li(p, pa)T 


2.821 


destroy ''Li 


'^Li{d, da)T 


3.175 


destroy ''Li 



Now we switch to examine the energy spectrum of BBNCRs. In principle 
the BBNCR energy spectrum will be determined by acceleration mechanism 
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and we are lack of knowledge on this. In the realistic case, the energy of BB- 
NCRs should exceed 1.8 MeV but the amount is required to decrease quickly, 
especially above 3.337 MeV, namely the energy threshold of D(p, n)2H. Mo- 
tivated by the power law of cosmic rays today, we assume the spectrum of 
BBNCRs obeys a power law with index 4 from 2 MeV to 4 MeV ^. Below 2 
MeV we choose the power index to be 2 or even in uniform distribution 
like white noise. Our choice of energy spectrum of BBNCRs is depicted in 
figure 1. 




Energy(MeV) 

Figure 1: BBNCR energy spectrum. 
a is the power index (see eq. (5) for detail). 

^ The energy range is so narrow that small variation of the power index will not change 
the results too much. It is obvious that the larger the power index is, the less important 
the higher energy particles are. 

^ Exactly it is the energy range between 0.09 MeV and 2 MeV. In fact, BBNCRs come 
from SBBN particles, and the lower limit of the energy range of BBNCRs is near or a 
little above the temperature of SBBN particles. In the SBBN, ^Be is produced mainly 
during the epoch when the temperature is between 0.08 MeV and 0.04 MeV, so it seems 
reasonable that we take 0.09 MeV as the lower limit. Note that ^Be abundance is not 
sensitive to the change of the lower limit. 
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BBNCRs can affect element abundance via endothermic processes. Be- 
sides these BBNCRs will also involve into exothermic SBBN reactions. The 
exothermic processes are thought much less important than endothermic ones 
in the SBBN relatively. However via exothermic processes the non-thermal 
BBNCRs may have important inffuence upon the element abundances, espe- 
cially through which D abundance may be decreased. Complete classification 
of exothermic processes is summarized in Appendix A. Reactions added in 
our numerical computation following: 

1. Endothermic reactions, including the most important process to de- 
stroy ^Be - ''Be(p,pQ;)^He and a weighty origin of ^Be - ^Li(p, n)^Be; 

2. Exothermic reactions, including the collisions between isotopes of hy- 
drogen, especially destroying D - D(p, 7)^He, D{d,p)T, D{d, n)^He, and 
T{d,n)'^B.e, and the ones producing or destroying ''Li by energetic hy- 
drogen among the most relevant reactions for BBN [45] - ^He(t, 7)'^Li, 
and ^Li(p, a)'^He; 

3. The ones to produce and destroy ^Li, including exothermic reactions 
^He((i, 7)^Li and ^Li(p, a)'^He and endothermic reactions ^Be((i,^ He)^Li 
and ^Li((i, t)^Li. Note that the computation for ^Li abundance is nei- 
ther complete nor accurate and we will come back to this point in 
section 4. 

The above reactions are listed in table 2: 
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Table 2: Additional reactions added in our element abundance computation. 



Process 


i nresnold / Me V 


Lross oection 


D(p,7)3He 





[46] 


D{d,p)T 
D(c^,n)3He 





[47] 





[48] 


T(d,n)^He 
^He(ci,7)6Li 





[49] 





[50] 


4He(t,7)^Li 





[51] 





[52, 53] 


^Li(p, a)^He 
^Li(j9, n)^Be 





[53, 54] 


1.880 


[55] 


'^Be(p, pa)^He 
7Li(d, t)6Li 
^Be(rf,3He)6Li 


1.814 


[56] 


1.278 


[57] 


0.144 


[58] 



3 Element abundances after including BBN- 
CRs contributions 

Now we turn to element abundance computation after including BBNCRs 
contributions. Besides the standard contributions in SBBN, extra contribu- 
tions come from the processes where the thermal particles collide with non- 
thermal BBNCRs particles. According to the Boltzmann equation, variation 
of the abundance of nuclide i is described as 

-^ = -Y.'^^M + J2'^kmi (1) 

j k,i 

where [ij] is the rate for destroying nuclide i and [kl] is the rate for creating 
i,Yi = Xi/ Ai'i^ the abundance of nuclide i with Xi the mass fraction and Ai 
the mass number of nuclide i. The sum over j goes through all reactions to 
destroy nuclide i and the sum over fc, / goes through all reactions to produce 
nuclide i. The rate [ij] is defined by 

M^NaP^. (2) 

where Na is the Avogadro's number, p is the baryon energy density, and 
(ij) = (cTf Here a is the cross section of the reaction {ij — )■ kl), and v is 
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the relative velocity between these two particles i and j. The () stands for 
the mean over different relative velocities [59]. In the SBBN, () means the 
thermal average, while in the case with BBNCRs it can be computed as 

{(Tv)i2{T,a) = — J dcos^^— X / fi{Ei,T)dEiX 



— oo 



1 



f2{E2,a)dE2a{Ei)v{EuE2,cose). (3) 

^2 Jom 

Here distribution of SBBN particles f\{Ei,T) is the normalized Boltzmann 
distribution, 

where k is the Boltzmann's constant, T is the Universe temperature, and Ei 
is the energy of the SBBN particle. Thus the normalization constant Ki = 1 
and the energy range is (— oo, +oo). Distribution of BBNCRs is power law 
with index a (see figure 1), 

/2(£;2,«)0CE2"", (5) 

f2 or for 0.09 MeV < E2 < 2 MeV 
a = < , (6) 

[4 for 2 MeV < ^2 < 4 MeV 

where E2 is the energy of the BBNCR particle, with a = 2 for power law 
case and a = for uniform distribution case. We normalize the function /2 
from 2 MeV to 4 MeV, 

= I f2iE2,a)dE2, (7) 

and make /2 continuous at 2 MeV point. When the mass of the SBBN 
particle is noted as mi and that of the BBNCR as m2, we can compute the 
relative velocity with the angle of incidence 6, 

V =\vi - V2\ 



I2E, , 2^2 ^ [E^E; ^ (8) 

1 4:\ cosfe'. 



mi m2 V rnim2 
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and incident energy Ei 

Ei = ^rriiV^, (9) 

where is the mass of the incident particle. The normahzation constant 
over 6 is = dcos^ = 2. The cross sections in eq. (3) are taken from 
experimental data source DataBase EXFOR [60] and ENDF[61]. As a result, 
the incident particle is not simply the energetic particle, but on the criterion 
from nuclear experiments. Taking ^He(t,7)^Li for example, data of cross 
sections are obtained under the situation where ^He is taken as the incident 
particle. 

We calculate abundances utilizing the updated version [62, 63] of the 
Wagoner code [64] from [65] with appropriate modification in order to include 
new contributions from BBNCRs. The baryon-to-photon ratio is updated to 
6.23x 10~^°[66]. The rates of endothermal and exothermic processes in table 2 
where BBNCRs particles involve are multiplied by e. The new contributions 
are added to the code as new channels. Details of the cross sections we 
adopted can be found in Appendix B. 

As the cross section of ^Be(p, pa)^He is unknown from experiments, we 
take D(p, n)2H with a shift of the threshold energy as a substitute. If the 
difference is considerable, e.g., the D(p, n)2H cross section is x times the 
'^Be{p,pa)^B.e cross section, we need to replace e as xe. 

Now we turn to show our numerical results for element abundance after 
including new contributions from BBNCRs. Figure 2 and 3 show processing 
rates, which is defined as eYiVj[ij]/ H [67] {H is the expansion rate), of produc- 
ing and destroying ''Li/ ''Be. From the figures we can clearly see contributions 
from non-thermal processes, especially in the low temperature regime. 

^Be abundances under a = 2 for 0.09 MeV < E2 < 2MeV and different 
values of e are shown in figure 4, among which the red solid line stands for 
SBBN result (e = 0). We can see that e = 7 x 10"^ can destroy 70% ^Be. 
As a price, 5% D is destroyed (as a comparison, the number shifts to 1% for 
uniform distribution). From the figure we can see clearly that new contribu- 
tions from BBNCRs can account for ^Li abundance quite satisfactory. Note 
that for the processes whose reactants are both hydrogen, the effect on D 
abundance need to be doubled according to the symmetry. 

In figure 5 and 6 shows abundances of ''Li/ ''Be and ^Li, as well as all 
elements as a function of temperature respectively. Here the solid lines stand 
for the SBBN results, and the dashed lines for results with BBNCRs under 
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Figure 2: Processing rate of producing ^Li/^Be as a function of temperature. 



a = 2 for 0.09 MeV < E2 < 2 MeV and e = 7 x 10-^. While the BBNCRs 
can account for ^Li abundance, for the ^Li abundance, energetic BBNCRs 
can enhance ^Li one order of magnitude, but not enough compared to ob- 
servations. From the curves we can also see that the shifts of other element 
abundances after including BBNCRs contributions are usually tiny. 

For completeness we also show processing rates of producing and destroy- 
ing ^Li in figure 7 and 8. It is not hard to see that non-thermal contributions 
are quite similar to the ^Li case. BBNCRs play important role mainly when 
the Universe temperature falls below 0.04 MeV. 
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Figure 3: Processing rate of destroying ''Li/ ''Be as a function of temperature. 
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Figure 4: ''Be abundance as a function of temperature. 
Here e is taken as 0, 5 x 10~^, 7 x 10~^ and 1 x 10~^ respectively. 
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Figure 5: Abundances of ^Li/^Be and ^Li 
function of temperature with e = 7 x 10 . 
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Figure 6: Abundances of all elements 
function of temperature with e = 7 x 10 
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Figure 7: Processing rate of producing ^Li as a function of temperature. 
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Figure 8: Processing rate of destroying ^Li as a function of temperature. 
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4 Conclusions and discussions 



In this paper, we investigated whether cosmic rays in the BBN epoch (BB- 
NCRs) can account for Lithium problem or not. In order to keep the success 
of SBBN, the flux, energy range and spectrum of BBNCRs are severely con- 
strained. In the allowed parameter space, extra contributions from BBNCRs 
to ^Li abundance can fill the discrepancy between SBBN prediction and ob- 
servations. However BBNCRs can lift ^Li abundance in the SBBN an order 
of magnitude, but still less than that of observations. 

We need to point out some factors beyond investigations in the pa- 
per. First, measurements of cross sections of nuclear reactions, especially 
''Be(p,pa)^He is critical. Lack of knowledge on T{t, 7)^He (^He will decay to 
^Li in 0.8s) and ^He(t, 7)^Li cross sections brings uncertainty on ^Li results, 
too. If their cross section between 2 MeV and 4 MeV is about 0(mb), the 
contributions to ^Li abundance is comparable to that of ^He((i, 7)^Li. 

Second, we consider hydrogen as the only component of BBNCRs in this 
paper, however the possible ^He or ^He as BBNCRs particles will change our 
results. We can imagine that ^He and ^He are accelerated to the same energy 
and spectrum (in fact their energy should be lower and the spectrum be softer 
than hydrogen). In this case, effects of some reactions (e.g. ^He((i, 7)^Li and 
^He(t,7)^Li) will be doubled. At the same time, new reactions are triggered 
by energetic '^He and ^He (see table 3). The reaction ^Be(Q;, p)^°B will destroy 
''Be and produce ^°B. 

In our analysis, it is important that the amount of BBNCRs should be 
very low above 3.337 MeV, namely the D(p, n)2H threshold. In order to 
test the effect of the spectrum choice on the deuteron abundance, we tried 
different spectrums, and show some representative ones in figure 9. For each 
spectrum, we give the necessary fraction of BBNCRs, namely e, to destroy 
70% ^Be, and the corresponding percentage of destructed deuterium. The 
results are shown in table 4. The upper limit of BBNCRs spectrum is taken 
to be 10 MeV We find that BBNCRs above 3.337 MeV would destroy 
deuterium inevitably. Therefore, the main part of the realistic spectrum 
should lie between 1.814 MeV and 3.337 MeV, and a rapid drop below 1.814 
MeV and above 3.337 MeV is favorable. 

At last but not least, the critic issue for the BBNCRs is about acceleration 
mechanism during BBN and whether the required flux and energy range as 

We ignore endothermic reactions whose threshold is higher 
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Figure 9: Energy spectrum test. 
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Table 3: List of new reactions triggered by energetic ^He and ^He. 



Process Threshold/Me V 

6Li(3He, 2p)^Li 0.703 

^Li(3He,t)^Be 1.259 

^Be(a,p)i°B 1.800 

6Li(3He, rfa)3He 2.213 

^Be(3He,a)23He 2.268 

6Li(a, c/)2^He 2.455 

^Be(a,2a)3He 2.491 

^U{a,dfBe 2.608 

^Be(a,3He)SBe 2.635 

6Li(3He, rf)^B 2.965 

6Li(3He,np)^Be 3.171 



Table 4: The effect of the spectrum choice on the deuteron abundance. 



Description of spectrum 


Curves in figure 9 


e fraction 


D destroyed 


uniform distribution between 2 MeV and 4 MeV 


black, solid 


2x10-5 


1.9% 


uniform distribution between 0.09 MeV and 2 MeV, 
E-'^ between 2 MeV and 10 MeV 


black, dashed 


3x10-5 


10.6% 


E^e'^^ between 0.09 MeV and 10 MeV 


red, solid 


3x10-"* 


15% 


E'^e-^^ between 0.09 MeV and 10 MeV 


orange, solid 


LSxlO--* 


1.8% 


£;4g-(B/i.8i4)a between 0.09 MeV and 10 MeV 


blue, dashed 


1x10-"* 


1.2% 



well as spectrum can be induced. Anyhow the energetic particles are thought 
to be easily thermalized by photons and nuclei around before they collide with 
other nuclei [41, 42]. 
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A Classification of Exothermic Reactions 

Exothermic reactions where BBNCRs involve are listed in two tables. 
Table 5 shows SBBN reactions which are included in the Wagoner code and 
table 6 shows reactions besides those. 
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Table 5: Exothermic reactions - Part I 



Process Effect 

H(n, 7)D produce D, but not important, for n is too little 

D(n, 7)T destroy D, but not important, for n is too little 

D(p, 7)^He destroy D, see table 2 

D((i, n)^He destroy D, see table 2 

D((i, p)T destroy D, see table 2 

T(p, 7)^He destroy T""' 

T((i, n)^He destroy D, see table 2; destroy T, see the footnote 

^He((i,p)^He destroy D, but less important than D(p, 7)^He; 

destroy ^He, see the footnote 

^He((i, 7)^Li produce ^Li, see table 2 

^He(t,7)^Li produce ^Li, see table 2 

^Li(p, 7)''Be destroy ^Li, but less important than ^Li(p, a)'^He, for the 

electromagnetic cross section is smaller than the strong one 

^Li(p, Q;)^He destroy ^Li, see table 2 

^Li(p, a)'^He destroy ^Li, see table 2 

^Li((i, ?T.)2^He destroy ^Li, but less important than ''Li(j9, a)^He 

^Be((i,p)2^He destroy '''Be, but less important 

^Be(p, 7)®B destroy ^Be, but less important, for the electromagnetic 
cross section 



Table 6: Exothermic reactions - Part II 



Process Effect 

D((i, 7)^He destroy D, but less important than D((i, p)T or 

D((i, n)'^He, for the electromagnetic cross section 

T(i:,2n)'^He destroy T, but change little compared to the SBBN situation 

T(t,7)^He produce ^Li, see Section 4 

^He(t, 7)^Li produce ^Li, see Section 4 

^He(t, (i)*^He produce D, but change little compared to the SBBN situation 

^He(t,np)^He destroy T and ^He, but change little compared to the SBBN situation 

To be continued. . . 



T will decay to "^He, and observations only set the upper limit of ■^He[68] 
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Process Effect 



^Li{d, a)'^He destroy ^Li, but less important 
^Li((i,p)^Li destroy ^Li, but less important 

^Li((i, ?7,)^Be destroy ^Li, but less important 
^Li((i,pt)^He destroy ^Li, but less important 
^Li((i, n^He)^He destroy ^Li, but less important 
^Li(t,7)^Be destroy ^Li, but less important; 

produce ^Be, but less important than ^Li(t, n)^Be, 

for the electromagnetic cross section 
^Li(t,n)2^He destroy ^Li, but less important 
^Li(t,(i)^Li destroy ^Li, but less important 

^Li(t,p)^Li destroy ^Li, but less important 

^Li((i, 7)^Be produce ^Be, but less important than ^Li(t, n)^Be, 

for the electromagnetic cross section 
^Li(t,7)^°Be produce i°B, (^^Be decays to ^°B) 

but less important than '^Be{a,py^B 
'^Li(t,n)^Be produce ^Be, maybe important, constraints from observations [69] 
^Li(t,a)^He produce ^Li, but less important than T(t,7)^He or ^He(t,7)^Li 
^Li(t, 2n)2^He destroy ^Li, but less important 
'^Be()f:, 7)^°B produce ^°B, but less important than ^Be(a,p)^°B 
^Be(t,a)^Li produce ^Li, but less important than T(t,7)^He or ^He(t,7)^Li 
^Be(i, (i)2^He destroy ''Be, but less important 

^Be(t,p)^Be produce ^Be, maybe important, constraints from observations [69] 
^Be(t, np)2'^He destroy ''Be, but less important 

^Be(t,^ He)^Li transformation between ^Li and ''Be, but less important than 
''Li(p, n) ''Be 



B Cross Sections and Details of Computation 
Scheme 

Cross sections of the reactions included in our computation are shown in 
figure 10-21. The possible errors here can be accommodated by changing the 
amount of BBNCRs, namely varying the free parameter e. 

Cross sections below the first available point from the experiment is set 
to zero, and above the last available point the cross sections are set to be 
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flat. Among the experimental points, cross sections are interpolated simply 
linearly. However for D(p, 7)'^He, '^Li{p,a)'^}ic, ^Li(p, a)^He and T((i, n)'^He, 
the cross sections are interpolated double-logarithmic linearly, in order to 
make them more smooth. 

The integral limit over thermal spectrum is not infinite. Instead we use 
O.OOIT as the lower hmit and 8T as the upper hmit, where T is the Universe 
temperature. Integration is carried out in six-point Gaussian scheme. 
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EHDF Keciuest 529, 2ail-Sep-3a, 89:49: 52 
EXFOP R^quEst: 1515/1, 2ail-Sep-3a a9:48:25 



10~^ 10~* 1 IQ IQ^ 

EHDF/B-UII.a: H-2(P,C)HE-3 




10~^ 10~* 1 IQ IQ^ 

Incident Energy (MeU) 



Figure 10: Cross sections of D(p, 7)^He 



l-H-2<D,P>l-H-3 
EXFaP Request: 1516/1, 2ail-Sep-3a a9:53:58 



10~^ 10~* 1 IQ 




10~^ 10~* 1 IQ 

Incident Energy (MeU) 



Figure 11: Cross sections of D{d,p)T 
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l-H-2(D,H>2-HE-3 
EXFOP R^quEst: 1517/1, 2ail-Sep-3a 09:56:22 




Figure 12: Cross sections of D((i, n)^He 




Figure 13: Cross sections of T{d,n)^}ie 
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l-H-2(A,G>3-LI-6 
EXFOP R^quEst: 1522/1, 2ail-Sep-3a 10:05:39 



-n- 1? 

" 'I 


91 Mohf X4 : Aa54BBB3 X^-B(us : 1, pt : 21, e 
Bl Rober-tson X4 1^1190002 /2=B(us : 1, pt : 14, e 


(■r':4-29K) 


■' 







ID 

Incident Energy (MeU) 



Figure 14: Cross sections of ^B.e{d, •y^Li 



l-H-3(A,G>3-LI-7 
EXFOR Request: 1523/1, 2011-Sep-30 10:07:35 

Q .B Q .8 1 .Q 1 .a 




□ .6 0.8 1 .□ 1.2 

Incident Energy (MeU) 



Figure 15: Cross sections of '^He(t, 7)^Li 
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3-L I -5 ( P , HE3 ) 2-HE-4 
EXFOP R^quEst: 1524/1, 2ail-Sep-30 19:10:19 



Int^r'Pola'tQd i 



-n- 1982 Szabo 
-O 1979 Eluyn 
1974 Gould 
-> 1967 Fiedler 
-if 1966 Ceneinhrii-dt X4 



:FBBeiBB3 ?:^-B(us:l,pt:5,dr'r:l5%) 

: FBB12BBB ?^^=B . 211656(us : iLpt : 33, err- : 9J!> 

:fll417BB3 /^=B.862189(UE Jsjpt:58,eri':l-12Si) 

: ni393BB2 jf*=2 . 72B19CuE : llfflt : 5, err- : 15-25!« ) 

:ni367002 X'=1.19265(us:l,^t:l5,eiT:4-lBX> 




Incident Energy (MeU) 



Figure 16: Cross sections of ^Li(p, a)^He 



3-LI-7<P,A)2-HE-4 
EXFOP Request: 1525/1, 2ail-Sep-3a 10:13:18 

10~^ 1 



Interpolated curue 

^ 1967 Fiedler X4:ftl393BB3 ?:^-B . 166544(us : 1, pt : 7, s 

-n- 2BB1 Spraher X4:CB619BB5 ?^^=B(u3:l,pt:l3,ei 

- 1986 Rolfs X4:Ca417BB2 ?^^=B<us : 1, pt : 5, er 

1964 Mani ]{4:ni465BB2 ;<'=B(us : 1, pt : 42j 



Incident Energy (MeU) 



Figure 17: Cross sections of ^Li(p, Q;)^He 
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3-LI-7<P,H)4-BE-7 
EXFOP R^quEst: 1532/1, 2ail-Sep-3a 10:23:15 



/^-1.695ZZ(u3:l,pt: 141, err-: 5-14%) 
?^^=2.a2277<us:l,pt:73,Gr'r:6K> 
?^^=12 . 5854(us : 1, pt : 11, Gcr : IX) 




Incident Energy (MeU) 



Figure 18: Cross sections of ^Li(p, n)^Be 



1-H-2(P,H+P)1-H-1 
EXFOP Request: 1534/1, 2ail-Sep-30 19:25:13 




Figure 19: Cross sections of D(p, r;,)2H 
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3-LI-7<D,T)3-LI-G 
EXFOP R^quEst: 1535/1, 2ail-Sep-3a 10:28:64 



2 3 4 




Incident Energy (MeU) 

Figure 20: Cross sections of ^Li(c?, t)^Li 



EHDF Keciuest 524, 2eil-Sep-3e,10 : 30: 25 

a 2Q 4Q BO 

• TEHDL-2099: BE-7CD,HE3)LI-6 

- Interpolated curve 




□ 20 40 BO 

Incident Energy (MeU) 



Figure 21: Cross sections of ^Be((i,'^ He)®Li 
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